\(\int \frac {\sec ^2(e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx\) [241]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [B] (warning: unable to verify)
   Fricas [A] (verification not implemented)
   Sympy [F]
   Maxima [F]
   Giac [F(-2)]
   Mupad [F(-1)]

Optimal result

Integrand size = 35, antiderivative size = 124 \[ \int \frac {\sec ^2(e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=-\frac {\sqrt {2} \arctan \left (\frac {\sqrt {a} \tan (e+f x)}{\sqrt {2} \sqrt {a+a \sec (e+f x)}}\right )}{\sqrt {a} (c-d) f}+\frac {2 c \arctan \left (\frac {\sqrt {a} \sqrt {d} \tan (e+f x)}{\sqrt {c+d} \sqrt {a+a \sec (e+f x)}}\right )}{\sqrt {a} (c-d) \sqrt {d} \sqrt {c+d} f} \]

[Out]

-arctan(1/2*a^(1/2)*tan(f*x+e)*2^(1/2)/(a+a*sec(f*x+e))^(1/2))*2^(1/2)/(c-d)/f/a^(1/2)+2*c*arctan(a^(1/2)*d^(1
/2)*tan(f*x+e)/(c+d)^(1/2)/(a+a*sec(f*x+e))^(1/2))/(c-d)/f/a^(1/2)/d^(1/2)/(c+d)^(1/2)

Rubi [A] (verified)

Time = 0.43 (sec) , antiderivative size = 124, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.143, Rules used = {4061, 3880, 209, 4052, 211} \[ \int \frac {\sec ^2(e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\frac {2 c \arctan \left (\frac {\sqrt {a} \sqrt {d} \tan (e+f x)}{\sqrt {c+d} \sqrt {a \sec (e+f x)+a}}\right )}{\sqrt {a} \sqrt {d} f (c-d) \sqrt {c+d}}-\frac {\sqrt {2} \arctan \left (\frac {\sqrt {a} \tan (e+f x)}{\sqrt {2} \sqrt {a \sec (e+f x)+a}}\right )}{\sqrt {a} f (c-d)} \]

[In]

Int[Sec[e + f*x]^2/(Sqrt[a + a*Sec[e + f*x]]*(c + d*Sec[e + f*x])),x]

[Out]

-((Sqrt[2]*ArcTan[(Sqrt[a]*Tan[e + f*x])/(Sqrt[2]*Sqrt[a + a*Sec[e + f*x]])])/(Sqrt[a]*(c - d)*f)) + (2*c*ArcT
an[(Sqrt[a]*Sqrt[d]*Tan[e + f*x])/(Sqrt[c + d]*Sqrt[a + a*Sec[e + f*x]])])/(Sqrt[a]*(c - d)*Sqrt[d]*Sqrt[c + d
]*f)

Rule 209

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[b, 2]))*ArcTan[Rt[b, 2]*(x/Rt[a, 2])], x] /;
 FreeQ[{a, b}, x] && PosQ[a/b] && (GtQ[a, 0] || GtQ[b, 0])

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 3880

Int[csc[(e_.) + (f_.)*(x_)]/Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)], x_Symbol] :> Dist[-2/f, Subst[Int[1/(2
*a + x^2), x], x, b*(Cot[e + f*x]/Sqrt[a + b*Csc[e + f*x]])], x] /; FreeQ[{a, b, e, f}, x] && EqQ[a^2 - b^2, 0
]

Rule 4052

Int[(csc[(e_.) + (f_.)*(x_)]*Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)])/(csc[(e_.) + (f_.)*(x_)]*(d_.) + (c_)
), x_Symbol] :> Dist[-2*(b/f), Subst[Int[1/(b*c + a*d + d*x^2), x], x, b*(Cot[e + f*x]/Sqrt[a + b*Csc[e + f*x]
])], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0]

Rule 4061

Int[csc[(e_.) + (f_.)*(x_)]^2/(Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)]*(csc[(e_.) + (f_.)*(x_)]*(d_.) + (c_
))), x_Symbol] :> Dist[-a/(b*c - a*d), Int[Csc[e + f*x]/Sqrt[a + b*Csc[e + f*x]], x], x] + Dist[c/(b*c - a*d),
 Int[Csc[e + f*x]*(Sqrt[a + b*Csc[e + f*x]]/(c + d*Csc[e + f*x])), x], x] /; FreeQ[{a, b, c, d, e, f}, x] && N
eQ[b*c - a*d, 0] && (EqQ[a^2 - b^2, 0] || EqQ[c^2 - d^2, 0])

Rubi steps \begin{align*} \text {integral}& = -\frac {\int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)}} \, dx}{c-d}+\frac {c \int \frac {\sec (e+f x) \sqrt {a+a \sec (e+f x)}}{c+d \sec (e+f x)} \, dx}{a (c-d)} \\ & = \frac {2 \text {Subst}\left (\int \frac {1}{2 a+x^2} \, dx,x,-\frac {a \tan (e+f x)}{\sqrt {a+a \sec (e+f x)}}\right )}{(c-d) f}-\frac {(2 c) \text {Subst}\left (\int \frac {1}{a c+a d+d x^2} \, dx,x,-\frac {a \tan (e+f x)}{\sqrt {a+a \sec (e+f x)}}\right )}{(c-d) f} \\ & = -\frac {\sqrt {2} \arctan \left (\frac {\sqrt {a} \tan (e+f x)}{\sqrt {2} \sqrt {a+a \sec (e+f x)}}\right )}{\sqrt {a} (c-d) f}+\frac {2 c \arctan \left (\frac {\sqrt {a} \sqrt {d} \tan (e+f x)}{\sqrt {c+d} \sqrt {a+a \sec (e+f x)}}\right )}{\sqrt {a} (c-d) \sqrt {d} \sqrt {c+d} f} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.72 (sec) , antiderivative size = 141, normalized size of antiderivative = 1.14 \[ \int \frac {\sec ^2(e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=-\frac {2 \left (\sqrt {d} \sqrt {c+d} \arctan \left (\frac {\sin \left (\frac {1}{2} (e+f x)\right )}{\sqrt {\cos (e+f x)}}\right )-\sqrt {2} c \arctan \left (\frac {\sqrt {2} \sqrt {d} \sin \left (\frac {1}{2} (e+f x)\right )}{\sqrt {c+d} \sqrt {\cos (e+f x)}}\right )\right ) \cos \left (\frac {1}{2} (e+f x)\right )}{(c-d) \sqrt {d} \sqrt {c+d} f \sqrt {\cos (e+f x)} \sqrt {a (1+\sec (e+f x))}} \]

[In]

Integrate[Sec[e + f*x]^2/(Sqrt[a + a*Sec[e + f*x]]*(c + d*Sec[e + f*x])),x]

[Out]

(-2*(Sqrt[d]*Sqrt[c + d]*ArcTan[Sin[(e + f*x)/2]/Sqrt[Cos[e + f*x]]] - Sqrt[2]*c*ArcTan[(Sqrt[2]*Sqrt[d]*Sin[(
e + f*x)/2])/(Sqrt[c + d]*Sqrt[Cos[e + f*x]])])*Cos[(e + f*x)/2])/((c - d)*Sqrt[d]*Sqrt[c + d]*f*Sqrt[Cos[e +
f*x]]*Sqrt[a*(1 + Sec[e + f*x])])

Maple [B] (warning: unable to verify)

Leaf count of result is larger than twice the leaf count of optimal. \(503\) vs. \(2(101)=202\).

Time = 18.95 (sec) , antiderivative size = 504, normalized size of antiderivative = 4.06

method result size
default \(-\frac {\sqrt {-\frac {2 a}{\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}}\, \sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\, \left (2 \sqrt {\left (c +d \right ) \left (c -d \right )}\, \ln \left (\csc \left (f x +e \right )-\cot \left (f x +e \right )+\sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\right ) \sqrt {\frac {d}{c -d}}-c \sqrt {2}\, \ln \left (-\frac {2 \left (\sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\, \sqrt {2}\, \sqrt {\frac {d}{c -d}}\, c -\sqrt {2}\, \sqrt {\frac {d}{c -d}}\, \sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\, d +\sqrt {\left (c +d \right ) \left (c -d \right )}\, \left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right )-c +d \right )}{-c \left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right )+\left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right ) d +\sqrt {\left (c +d \right ) \left (c -d \right )}}\right )+c \sqrt {2}\, \ln \left (-\frac {2 \left (-\sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\, \sqrt {2}\, \sqrt {\frac {d}{c -d}}\, c +\sqrt {2}\, \sqrt {\frac {d}{c -d}}\, \sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\, d +\sqrt {\left (c +d \right ) \left (c -d \right )}\, \left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right )+c -d \right )}{c \left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right )-\left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right ) d +\sqrt {\left (c +d \right ) \left (c -d \right )}}\right )\right )}{2 f a \sqrt {\left (c +d \right ) \left (c -d \right )}\, \left (c -d \right ) \sqrt {\frac {d}{c -d}}}\) \(504\)

[In]

int(sec(f*x+e)^2/(c+d*sec(f*x+e))/(a+a*sec(f*x+e))^(1/2),x,method=_RETURNVERBOSE)

[Out]

-1/2/f/a/((c+d)*(c-d))^(1/2)/(c-d)/(d/(c-d))^(1/2)*(-2*a/((1-cos(f*x+e))^2*csc(f*x+e)^2-1))^(1/2)*((1-cos(f*x+
e))^2*csc(f*x+e)^2-1)^(1/2)*(2*((c+d)*(c-d))^(1/2)*ln(csc(f*x+e)-cot(f*x+e)+((1-cos(f*x+e))^2*csc(f*x+e)^2-1)^
(1/2))*(d/(c-d))^(1/2)-c*2^(1/2)*ln(-2*(((1-cos(f*x+e))^2*csc(f*x+e)^2-1)^(1/2)*2^(1/2)*(d/(c-d))^(1/2)*c-2^(1
/2)*(d/(c-d))^(1/2)*((1-cos(f*x+e))^2*csc(f*x+e)^2-1)^(1/2)*d+((c+d)*(c-d))^(1/2)*(-cot(f*x+e)+csc(f*x+e))-c+d
)/(-c*(-cot(f*x+e)+csc(f*x+e))+(-cot(f*x+e)+csc(f*x+e))*d+((c+d)*(c-d))^(1/2)))+c*2^(1/2)*ln(-2*(-((1-cos(f*x+
e))^2*csc(f*x+e)^2-1)^(1/2)*2^(1/2)*(d/(c-d))^(1/2)*c+2^(1/2)*(d/(c-d))^(1/2)*((1-cos(f*x+e))^2*csc(f*x+e)^2-1
)^(1/2)*d+((c+d)*(c-d))^(1/2)*(-cot(f*x+e)+csc(f*x+e))+c-d)/(c*(-cot(f*x+e)+csc(f*x+e))-(-cot(f*x+e)+csc(f*x+e
))*d+((c+d)*(c-d))^(1/2))))

Fricas [A] (verification not implemented)

none

Time = 0.57 (sec) , antiderivative size = 1041, normalized size of antiderivative = 8.40 \[ \int \frac {\sec ^2(e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\text {Too large to display} \]

[In]

integrate(sec(f*x+e)^2/(c+d*sec(f*x+e))/(a+a*sec(f*x+e))^(1/2),x, algorithm="fricas")

[Out]

[-1/2*(sqrt(2)*(a*c*d + a*d^2)*sqrt(-1/a)*log(-(2*sqrt(2)*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*sqrt(-1/a)*c
os(f*x + e)*sin(f*x + e) - 3*cos(f*x + e)^2 - 2*cos(f*x + e) + 1)/(cos(f*x + e)^2 + 2*cos(f*x + e) + 1)) - sqr
t(-a*c*d - a*d^2)*c*log(-((a*c^2 + 8*a*c*d + 8*a*d^2)*cos(f*x + e)^3 + a*d^2 + (a*c^2 + 2*a*c*d)*cos(f*x + e)^
2 - 4*sqrt(-a*c*d - a*d^2)*((c + 2*d)*cos(f*x + e)^2 - d*cos(f*x + e))*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))
*sin(f*x + e) - (6*a*c*d + 7*a*d^2)*cos(f*x + e))/(c^2*cos(f*x + e)^3 + (c^2 + 2*c*d)*cos(f*x + e)^2 + d^2 + (
2*c*d + d^2)*cos(f*x + e))))/((a*c^2*d - a*d^3)*f), -1/2*(sqrt(2)*(a*c*d + a*d^2)*sqrt(-1/a)*log(-(2*sqrt(2)*s
qrt((a*cos(f*x + e) + a)/cos(f*x + e))*sqrt(-1/a)*cos(f*x + e)*sin(f*x + e) - 3*cos(f*x + e)^2 - 2*cos(f*x + e
) + 1)/(cos(f*x + e)^2 + 2*cos(f*x + e) + 1)) - 2*sqrt(a*c*d + a*d^2)*c*arctan(2*sqrt(a*c*d + a*d^2)*sqrt((a*c
os(f*x + e) + a)/cos(f*x + e))*cos(f*x + e)*sin(f*x + e)/((a*c + 2*a*d)*cos(f*x + e)^2 - a*d + (a*c + a*d)*cos
(f*x + e))))/((a*c^2*d - a*d^3)*f), 1/2*(sqrt(-a*c*d - a*d^2)*c*log(-((a*c^2 + 8*a*c*d + 8*a*d^2)*cos(f*x + e)
^3 + a*d^2 + (a*c^2 + 2*a*c*d)*cos(f*x + e)^2 - 4*sqrt(-a*c*d - a*d^2)*((c + 2*d)*cos(f*x + e)^2 - d*cos(f*x +
 e))*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*sin(f*x + e) - (6*a*c*d + 7*a*d^2)*cos(f*x + e))/(c^2*cos(f*x + e
)^3 + (c^2 + 2*c*d)*cos(f*x + e)^2 + d^2 + (2*c*d + d^2)*cos(f*x + e))) + 2*sqrt(2)*(a*c*d + a*d^2)*arctan(sqr
t(2)*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*cos(f*x + e)/(sqrt(a)*sin(f*x + e)))/sqrt(a))/((a*c^2*d - a*d^3)*
f), (sqrt(a*c*d + a*d^2)*c*arctan(2*sqrt(a*c*d + a*d^2)*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*cos(f*x + e)*s
in(f*x + e)/((a*c + 2*a*d)*cos(f*x + e)^2 - a*d + (a*c + a*d)*cos(f*x + e))) + sqrt(2)*(a*c*d + a*d^2)*arctan(
sqrt(2)*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*cos(f*x + e)/(sqrt(a)*sin(f*x + e)))/sqrt(a))/((a*c^2*d - a*d^
3)*f)]

Sympy [F]

\[ \int \frac {\sec ^2(e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\int \frac {\sec ^{2}{\left (e + f x \right )}}{\sqrt {a \left (\sec {\left (e + f x \right )} + 1\right )} \left (c + d \sec {\left (e + f x \right )}\right )}\, dx \]

[In]

integrate(sec(f*x+e)**2/(c+d*sec(f*x+e))/(a+a*sec(f*x+e))**(1/2),x)

[Out]

Integral(sec(e + f*x)**2/(sqrt(a*(sec(e + f*x) + 1))*(c + d*sec(e + f*x))), x)

Maxima [F]

\[ \int \frac {\sec ^2(e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\int { \frac {\sec \left (f x + e\right )^{2}}{\sqrt {a \sec \left (f x + e\right ) + a} {\left (d \sec \left (f x + e\right ) + c\right )}} \,d x } \]

[In]

integrate(sec(f*x+e)^2/(c+d*sec(f*x+e))/(a+a*sec(f*x+e))^(1/2),x, algorithm="maxima")

[Out]

integrate(sec(f*x + e)^2/(sqrt(a*sec(f*x + e) + a)*(d*sec(f*x + e) + c)), x)

Giac [F(-2)]

Exception generated. \[ \int \frac {\sec ^2(e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\text {Exception raised: TypeError} \]

[In]

integrate(sec(f*x+e)^2/(c+d*sec(f*x+e))/(a+a*sec(f*x+e))^(1/2),x, algorithm="giac")

[Out]

Exception raised: TypeError >> an error occurred running a Giac command:INPUT:sage2:=int(sage0,sageVARx):;OUTP
UT:index.cc index_m i_lex_is_greater Error: Bad Argument Value

Mupad [F(-1)]

Timed out. \[ \int \frac {\sec ^2(e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\int \frac {1}{{\cos \left (e+f\,x\right )}^2\,\sqrt {a+\frac {a}{\cos \left (e+f\,x\right )}}\,\left (c+\frac {d}{\cos \left (e+f\,x\right )}\right )} \,d x \]

[In]

int(1/(cos(e + f*x)^2*(a + a/cos(e + f*x))^(1/2)*(c + d/cos(e + f*x))),x)

[Out]

int(1/(cos(e + f*x)^2*(a + a/cos(e + f*x))^(1/2)*(c + d/cos(e + f*x))), x)